Computational Exploration of Potential CFTR Binding Sites for Type I Corrector Drugs.

Cystic fibrosis (CF) is a recessive genetic disease that is caused by mutations in the cystic fibrosis transmembrane conductance regulator (CFTR) protein. The recent development of a class of drugs called "correctors", which repair the structure and function of mutant CFTR, has greatly enhanced the life expectancy of CF patients. These correctors target the most common disease causing CFTR mutant F508del and are exemplified by the FDA-approved VX-809. While one binding site of VX-809 to CFTR was recently elucidated by cryo-electron microscopy, four additional binding sites have been proposed in the literature and it has been theorized that VX-809 and structurally similar correctors may engage multiple CFTR binding sites. To explore these five binding sites, ensemble docking was performed on wild-type CFTR and the F508del mutant using a large library of structurally similar corrector drugs, including VX-809 (lumacaftor), VX-661 (tezacaftor), ABBV-2222 (galicaftor), and a host of other structurally related molecules. For wild-type CFTR, we find that only one site, located in membrane spanning domain 1 (MSD1), binds favorably to our ligand library. While this MSD1 site also binds our ligand library for F508del-CFTR, the F508del mutation also opens a binding site in nucleotide binding domain 1 (NBD1), which enables strong binding of our ligand library to this site. This NBD1 site in F508del-CFTR exhibits the strongest overall binding affinity for our library of corrector drugs. This data may serve to better understand the structural changes induced by mutation of CFTR and how correctors bind to the protein. Additionally, it may aid in the design of new, more effective CFTR corrector drugs.

Antitumor T-cell Immunity Contributes to Pancreatic Cancer Immune Resistance.

Pancreatic ductal adenocarcinoma (PDAC) is the third leading cause of cancer death in the United States. Pancreatic tumors are minimally infiltrated by T cells and are largely refractory to immunotherapy. Accordingly, the role of T-cell immunity in pancreatic cancer has been somewhat overlooked. Here, we hypothesized that immune resistance in pancreatic cancer was induced in response to antitumor T-cell immune responses and that understanding how pancreatic tumors respond to immune attack may facilitate the development of more effective therapeutic strategies. We now provide evidence that T-cell-dependent host immune responses induce a PDAC-derived myeloid mimicry phenomenon and stimulate immune resistance. Three KPC mouse models of pancreatic cancer were used: the mT3-2D (Kras+/LSL-G12D; Trp53+/LSL-R172H; Pdx1-Cre) subcutaneous and orthotopic models, as well as the KP1 (p48-CRE/LSL-Kras/Trp53 flox/flox ) subcutaneous model. KPC cancer cells were grown in immunocompetent and immunodeficient C57BL/6 mice and analyzed to determine the impact of adaptive immunity on malignant epithelial cells, as well as on whole tumors. We found that induced T-cell antitumor immunity, via signal transducer and activator of transcription 1 (STAT1), stimulated malignant epithelial pancreatic cells to induce the expression of genes typically expressed by myeloid cells and altered intratumoral immunosuppressive myeloid cell profiles. Targeting the Janus Kinase (JAK)/STAT signaling pathway using the FDA-approved drug ruxolitinib overcame these tumor-protective responses and improved anti-PD-1 therapeutic efficacy. These findings provide future directions for treatments that specifically disable this mechanism of resistance in PDAC.

Carotenoids promote lateral packing and condensation of lipid membranes.

Carotenoids are pigment molecules that protect biomembranes against degradation and may be involved in the formation of functional bacterial membrane microdomains. Little is known on whether different types of carotenoids have different effects on the membrane or if there is any concentration dependence of these effects. In this work, we present results from molecular dynamics simulations of phospholipid bilayers containing different amounts of either ß-carotene or zeaxanthin. Both ß-carotene and zeaxanthin show the ability to laterally condense the membrane lipids and reduce their inter-leaflet interactions. With increasing concentrations, both carotenoids increase the bilayer thickness and rigidity. The results reveal that carotenoids have similar effects to cholesterol on regulating the behavior of fluid-phase membranes, suggesting that they could function as sterol substitutes and confirming their potential role in the formation of functional membrane domains.

CAMERRA: An analysis tool for the computation of conformational dynamics by evaluating residue-residue associations.

A computational method which extracts the dominant motions from an ensemble of biomolecular conformations via a correlation analysis of residue-residue contacts is presented. The algorithm first renders the structural information into contact matrices, then constructs the collective modes based on the correlated dynamics of a selected set of dynamic contacts. Associated programs can bridge the results for further visualization using graphics software. The aim of this method is to provide an analysis of conformations of biopolymers from the contact viewpoint. It may assist a systematical uncovering of conformational switching mechanisms existing in proteins and biopolymer systems in general by statistical analysis of simulation snapshots. In contrast to conventional correlation analyses of Cartesian coordinates (such as distance covariance analysis and Cartesian principal component analysis), this program also provides an alternative way to locate essential collective motions in general. Herein, we detail the algorithm in a stepwise manner and comment on the importance of the method as applied to decoding allosteric mechanisms. © 2018 Wiley Periodicals, Inc.

Characterizing protein conformations by correlation analysis of coarse-grained contact matrices.

We have developed a method to capture the essential conformational dynamics of folded biopolymers using statistical analysis of coarse-grained segment-segment contacts. Previously, the residue-residue contact analysis of simulation trajectories was successfully applied to the detection of conformational switching motions in biomolecular complexes. However, the application to large protein systems (larger than 1000 amino acid residues) is challenging using the description of residue contacts. Also, the residue-based method cannot be used to compare proteins with different sequences. To expand the scope of the method, we have tested several coarse-graining schemes that group a collection of consecutive residues into a segment. The definition of these segments may be derived from structural and sequence information, while the interaction strength of the coarse-grained segment-segment contacts is a function of the residue-residue contacts. We then perform covariance calculations on these coarse-grained contact matrices. We monitored how well the principal components of the contact matrices is preserved using various rendering functions. The new method was demonstrated to assist the reduction of the degrees of freedom for describing the conformation space, and it potentially allows for the analysis of a system that is approximately tenfold larger compared with the corresponding residue contact-based method. This method can also render a family of similar proteins into the same conformational space, and thus can be used to compare the structures of proteins with different sequences.

Effects of carotenoids on lipid bilayers.

Carotenoids have been found to be important in improving the integrity of biomembranes in eukaryotes. However, the molecular details of how carotenoids modulate the physical properties of biomembranes are unknown. To this end, we have conducted a series of molecular dynamics simulations of different biologically-relevant membranes in the presence of carotenoids. The carotenoid effect on the membrane was found to be specific to the identity of the carotenoid and the composition of the membrane itself. Therefore, different classes of carotenoids produce a different effect on the membrane, and different membrane phases are affected differently by carotenoids. It is apparent from our data that carotenoids do trigger the bilayer to become thinner. The mechanism by which this occurs depends on two competing factors, the ability of the lipid tails of opposing monolayers to either (1) compress or (2) interdigitate as the bilayer condenses. Indeed, carotenoids directly influence the physical properties via these two mechanisms, thus compacting the bilayer. However, the degree to which these competing mechanisms are utilized depends on the bilayer phase and the carotenoid identity.

DMSO enhanced conformational switch of an interfacial enzyme.

Interfacial proteins function in unique heterogeneous solvent environments, such as water-oil interfaces. One important example is microbial lipase, which is activated in an oil-water emulsion phase and has many important enzymatic functions. A unique aprotic dipolar organic solvent, dimethyl sulfoxide (DMSO), has been shown to increase the activity of lipases, but the mechanism behind this enhancement is still unknown. Here, all-atom molecular dynamics simulations of lipase in a binary solution were performed to examine the effects of DMSO on the dynamics of the gating mechanism. The amphiphilic α5 region of the lipase was a focal point for the analysis, since the structural ordering of α5 has been shown to be important for gating under other perturbations. Compared to the closed-gorge ensemble in an aqueous environment, the conformational ensemble shifts towards open-gorge structures in the presence of DMSO solvents. Increased width of the access channel is particularly prevalent in 45% and 60% DMSO concentrations (w/w). As the amount of DMSO increases, the α5 region of the lipase becomes more α-helical, as we previously observed in studies that address water-oil interfacial and high pressure activation. We believe that the structural ordering of α5 plays an essential role on gating and lipase activity.

The Promiscuity of Allosteric Regulation of Nuclear Receptors by Retinoid X Receptor.

The promiscuous protein retinoid X receptor (RXR) displays essential allosteric regulation of several members in the nuclear hormone receptor superfamily via heterodimerization and (anti)cooperative binding of cognate ligands. Here, the structural basis of the positive allostery of RXR and constitutive androstane receptor (CAR) is revealed. In contrast, a similar computational approach had previously revealed the mechanism for negative allostery in the complex of RXR and thyroid receptor (TR). By comparing the positive and negative allostery of RXR complexed with CAR and TR respectively, we reported the promiscuous allosteric control involving RXR. We characterize the allosteric mechanism by expressing the correlated dynamics of selected residue-residue contacts which was extracted from atomistic molecular dynamics simulation and statistical analysis. While the same set of residues in the binding pocket of RXR may initiate the residue-residue interaction network, RXR uses largely different sets of contacts (only about one-third identical) and allosteric modes to regulate TR and CAR. The promiscuity of RXR control may originate from multiple factors, including (1) the frustrated fit of cognate ligand 9c to the RXR binding pocket and (2) the different ligand-binding features of TR (loose) versus CAR (tight) to their corresponding cognate ligands.

Pressure-induced conformational switch of an interfacial protein.

A special class of proteins adopts an inactive conformation in aqueous solution and activates at an interface (such as the surface of lipid droplet) by switching their conformations. Lipase, an essential enzyme for breaking down lipids, serves as a model system for studying such interfacial proteins. The underlying conformational switch of lipase induced by solvent condition is achieved through changing the status of the gated substrate-access channel. Interestingly, a lipase was also reported to exhibit pressure activation, which indicates it is drastically active at high hydrostatic pressure. To unravel the molecular mechanism of this unusual phenomenon, we examined the structural changes induced by high hydrostatic pressures (up to 1500 MPa) using molecular dynamics simulations. By monitoring the width of the access channel, we found that the protein undergoes a conformational transition and opens the access channel at high pressures (>100 MPa). Particularly, a disordered amphiphilic α5 region of the protein becomes ordered at high pressure. This positive correlation between the channel opening and α5 ordering is consistent with the early findings of the gating motion in the presence of a water-oil interface. Statistical analysis of the ensemble of conformations also reveals the essential collective motions of the protein and how these motions contribute to gating. Arguments are presented as to why heightened sensitivity to high-pressure perturbation can be a general feature of switchable interfacial proteins. Further mutations are also suggested to validate our observations. Proteins 2016; 84:820-827. © 2016 Wiley Periodicals, Inc.

Investigation of Carbohydrate Recognition via Computer Simulation.

Carbohydrate recognition by proteins, such as lectins and other (bio)molecules, can be essential for many biological functions. Recently, interest has arisen due to potential protein and drug design and future bioengineering applications. A quantitative measurement of carbohydrate-protein interaction is thus important for the full characterization of sugar recognition. We focus on the aspect of utilizing computer simulations and biophysical models to evaluate the strength and specificity of carbohydrate recognition in this review. With increasing computational resources, better algorithms and refined modeling parameters, using state-of-the-art supercomputers to calculate the strength of the interaction between molecules has become increasingly mainstream. We review the current state of this technique and its successful applications for studying protein-sugar interactions in recent years.

Mapping allostery through computational glycine scanning and correlation analysis of residue-residue contacts.

Understanding allosteric mechanisms is essential for the physical control of molecular switches and downstream cellular responses. However, it is difficult to decode essential allosteric motions in a high-throughput scheme. A general two-pronged approach to performing automatic data reduction of simulation trajectories is presented here. The first step involves coarse-graining and identifying the most dynamic residue-residue contacts. The second step is performing principal component analysis of these contacts and extracting the large-scale collective motions expressed via these residue-residue contacts. We demonstrated the method using a protein complex of nuclear receptors. Using atomistic modeling and simulation, we examined the protein complex and a set of 18 glycine point mutations of residues that constitute the binding pocket of the ligand effector. The important motions that are responsible for the allostery are reported. In contrast to conventional induced-fit and lock-and-key binding mechanisms, a novel "frustrated-fit" binding mechanism of RXR for allosteric control was revealed.

Effects of branched O-glycosylation on a semiflexible peptide linker.

Glycosylation is an essential modification of proteins and lipids by the addition of carbohydrate residues. These attached carbohydrates range from single monomers to elaborate branched glycans. Here, we examine how the level of glycosylation affects the conformation of a semiflexible peptide linker using the example of the hinge peptide from immunoglobulin A. Three sets of atomistic models of this hinge peptide with varying degrees of glycosylation are constructed to probe how glycosylation affects the physical properties of the linker. We found that glycosylation greatly altered the predominant conformations of the peptide, causing it to become elongated in reference to the unglycosylated form. Furthermore, glycosylation restricts the conformational exploration of the peptide. At the residue level, glycans are found to introduce a bias for the formation of more extended secondary structural elements for glycosylated serines. Additionally, the flexibility of this semiflexible proline-rich peptide is significantly reduced by glycosylation.

Solvent-induced α- to 3(10)-helix transition of an amphiphilic peptide.

The amphiphilic peptide of the triacylglycerol lipase derived from Pseudomonas aeruginosa plays a critical role in guarding the gate for ligand access. Conformations of this peptide at several water-oil interfaces and in protein environments were compared using atomistic simulations with explicit solvents. In oil-containing solvents, this peptide is able to retain a folded structure. Interestingly, when the peptide is immersed in a low-polarity solvent environment, it exhibits a "coalesced" helix structure, which has both α- and 3(10)-helix components. The observation that the 3(10)-helical conformation is populated in a highly nonpolar environment is consistent with a previous report on polymethylalanine. Frequent interconversions of the secondary structure (between α-helix and 3(10)-helix) of the peptide are also observed. We further studied how this solvent-induced structural transition may be connected to the trigger mechanism of lipase gating and how the lipase senses the hydrophobic-hydrophilic interface.

Solvent-dependent gating motions of an extremophilic lipase from Pseudomonas aeruginosa.

Understanding how organic solvent-stable proteins can function in anhydrous and often complex solutions is essential for the study of the interaction of protein and molecular immiscible interfaces and the design of efficient industrial enzymes in nonaqueous solvents. Using an extremophilic lipase from Pseudomonas aeruginosa as an example, we investigated the conformational dynamics of an organic solvent-tolerant enzyme in complex solvent milieux. Four 100-ns molecular dynamics simulations of the lipase were performed in solvent systems: water, hexane, and two mixtures of hexane and water, 5% and 95% (w/w) hexane. Our results show a solvent-dependent structural change of the protein, especially in the region that regulates the admission of the substrate. We observed that the lipase is much less flexible in hexane than in aqueous solution or at the immiscible interface. Quantified by the size of the accessible channel, the lipase in water has a closed-gate conformation and no access to the active site, while in the hexane-containing systems, the lipase is at various degrees of open-gate state, with the immiscible interface setup being in the widely open conformation ensembles. The composition of explicit solvents in the access channel showed a significant influence on the conformational dynamics of the protein. Interestingly, the slowest step (bottleneck) of the hexane-induced conformational switch seems to be correlated with the slow dehydration dynamics of the channel.